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ABSTRACT 



Context. Recent observations with the spectropolarimeter on-board the HINODE satelhte have found abundant horizontal magnetic 
fields in the internetwork quiet Sun. 

Aims. We compare the results on the horizontal fields obtained at 630 nm with ground-based observations at 1.56 yum where the 
sensitivity to magnetic fields is larger than in the visible. 

Methods. We obtained 30 sec-integrated spectropolarimetric data of quiet Sun on disc centre during a period of extremely stable and 
good seeing. The data have a rms noise in polarization of around 2 • 10~^ of the continuum intensity. The low noise level allowed for 
an inversion of the spectra with the SIR code. We compare the inversion results with proxies for the determination of magnetic flux. 
Results. We confirm the presence of the horizontal fields in the quiet Sun internetwork as reported from the satellite data, with voids 
of some granules extent of nearly zero linear polarization signal. Voids in the circular polarization signal are only of granular scale. 
More than 60 % of the surface show polarization signals above four times the rms noise level. We find that the total magnetic flux 
contained in the more inclined to horizontal fields (y > 45°) is smaller by a factor of around 2 than that of the less inclined fields. The 
proxies for flux determination are seen to suff'er from a strong influence of the thermodynamic state of the atmosphere, and hence, 
seem to be unreliable. 

Conclusions. During spells of good seeing conditions, adaptive optics can render ground-based slit- spectrograph observations at a 
70-cm telescope equivalent to the seeing-free space-based data of half-meter class telescopes. We suggest that the diff'erence of the 
ratio of horizontal to transversal flux between the ground-based infrared data and the satellite-based visible data is due to the diff'erent 
formation heights of the respective spectral lines. We caution that the true amount of magnetic flux cannot be derived directly from 
the spectra. For purely horizontal flux, one would need its vertical extension that has to estimated by an explicit modeling with the 
observed spectra as boundary conditions, or has to be taken from MHD simulations. Time-series of the evolution of the magnetic flux 
and chromospheric diagnostics are needed to address its possible contribution to chromospheric heating. 

Key words. Sun: magnetic fields. Sun: photosphere 



1. Introduction 

Recent studies using observations of magnetic fields in quiet 
Sun regions obtai ned with the spectro polarimeter on-board the 
I HINODE satellite (iKosugi et al.l2007h have found abundant hor- 
. izontal magnetic fields torozco Suarez et alJ 120071: iLites et all 
■ |2QQ8[ ; OR07 and LI08 in the following), whose existence 
was predicted before from the mismatch of observed average 
magnetic flux between Ze eman and Hanle measurements by 
' iTrujillo Bueno et aP (l2004l) . Studies of photospheric magnetic 
fields using Zeeman- sensitive spectral lines in the visible (VIS) 
wavelength range, however, suffer from a severe drawback 
inherent to the wavelength: the thermal broadening (cx A) is 
dominating over the Zeeman splitting due to magnetic field 
(oc A^). For spectral lines in the near-infrared (IR), the relation 
is more favorable. The retrieval of magnetic field properties 
from spectra in the VIS thus is then less reliable, because the 
lines are in the weak-field limit that imposes severe limits 
(iMartinez Gonzalez et ani2006bl l2QQ8cl) . The separation of the 
oppositely polarized spectral components is not proportional 
to the field strength, and the information on field strength and 



magnetic flux both influence the amplitude of the polarization 
signal in similar ways. Observations of solar magnetic fields 
commonly do not resolve the magnetic structures, thus, the 
polarization signal amplitude, which is always measured 
relative to intensity, gets further complicated by the question 
of spatial resolution, stray light a nd the thermo dynamics of 
the solar atmosphere (Martin ez Gonzalez et alJ >2006b). An 
accurate determination of magnetic field properties from the 
weak polarization signals prevalent in the quiet Sun in VIS lines 
thus gets strongly dependent on the quality of the observations, 
on the signal-to-n oise ratio (Bellot Rubio & Collad osI 120031 
Rezaei et al]|2007l) and the spatial resolution (.Khomenko et alJ 



2005bl) . IR spectra at 1.56 fim allow to disentangle most of 
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the ambiguities because the field strength can be read oflT from 
the splitting of the polarization, or in some cases, even the 
intensity components of the spectral line for ni agnetic fields 
above around 400 G (.Beckll2006l:lBeck et al.ll2007b . It is thus not 
surprising that studies befo r e HIN O DE using solely VIS lines 
fean chez Alm eida & Lites 2000"; Dominguez Cerdena et alJ 
200 3]; jSocas-Na varro & Sanchez Almeida 120031 

Lite s & Socas-Navarrol 
l200j or IR lines 



iLin 



20041: 

Hall; 



J 



Socas-Navarro & Lite: 
Khomenko et al.i i2003l: 




Fig. 1. Overview maps of the observations on 21.5.08. Top: TIP Op. 001, taken at around UT 8:30. Bottom: TIP Op. 005, from 
around UT 9:30. Left to right: G-band image, 7^, L^ax, ^max^ polarity. Scale is in arcsec, grid lines have a 10 arcsec spacing. The 
blue diamonds denote the location of the profile shown in Fig.lH the red diamonds the locations of those shown in Appendix lAl 



'Martinez G onzalez et al.l l2006al: iDommguez C erdenaet 
t2006a) did not agree especially on the question of the intrinsic 
field strength: the VIS observations indicated a significant 
fraction of strong kG fields, whereas the IR was dominated by 
weaker hG fields. 

The contradictory findings from VIS or IR s pectral lines 
could only recently be reconciled to some degree. iRezaei et all 
(l2007l) retrieved for the first time a predominance of weak fields 
in the internetwork (IN) from observations at 630 nm in a study 
on the relation between photo spheric fields and chromo spheric 
emission. They used ground-based observations that were ob- 
tained using the Kiepenheuer-Institute adaptive optics system 
(von der Liihe et al. 2003) and the POLIS spectropolarimeter 
( Beck et al.i 2005b) at the German Vacuum Tower Telescope 
(VTT) in Izana. Martine z Gonzal ez et al.i (i2008c) showed that 
for simultaneous observations in VIS and IR lines taken under 
identical seeing conditions, the VIS observations are compatible 
with the weak fields retrieved from the corresponding IR spectra. 
The VIS spectra alone were, however, found again to be biased 
towards kG fields due to the limitations imposed by the weak- 
field limit. 

The new findings on the magnetic fields in the internetwork 
quiet Sun from the HINODE in VIS lines thus should again be 
best compared and confirmed by corresponding observations in 
the IR. During a recent observation campaign, we were able to 
obtain IR spectra at 1.56yum during a period of extremely good 
and stable seeing, yielding a nearly uniform spatial resolution in 
time and space that was near the diff'raction limit of the VTT at 
these wavelengths (0'.'6). We describe the observations in Sect.[2l 



and explain additional analysis steps and the inv ersion of the 
spectra with the SIR code ( Rui z Cobo & d el Toro Iniesta 199^ 
in Sect. [3l The results are presented in Sect. IH and summarized 
and discussed in Sect. \5\ Our conclusions are given in Sect. [6l 

2, Observations 

On 21/05/2008, we observed a region of quiet Sun on 
disc centre with the Tenerife Infrared Pol arimeter (TIP; 
iMartmez Pillet et al.|[T999l: ICollados et al.ll2007h at the German 
Vacuum Tower Telescope (VTT) in Izana. We used the 1.56 jim 
range that includes two Zeeman sensitive Fe I solar spectral lines 
at 1564.8 nm (geff = 3) and 1565.2 nm (ge/f = 1.5). The slit 
width of 100 yum corresponded to 0^.^45. The slit was stepped in 
0'.'5 steps across the solar image in the focal plane. The spatial 
sampling along the slit was 0'.'175; the spectra were binned later 
to a 0^.^35 sampling that is sufficient to sample the diff'raction 
limit of the VTT (0'.'6 at 1.5 jim). Due to the extremely good 
and especially stable seeing conditions, we used an integration 
time of 30 sees per scan step. Two maps of 40''x70'' were taken 
from UT 7:58 until 8:43 (Operation 001), and from UT 9:16 un- 
til 10:01 (Operation 005 ). The TIP d ata were reduced with the 
respective routines (e.g., ISchlichenmaier & Collado s 2002) in- 
cluding the correction for the telescope polarization dBeck et al.l 
2005a). The noise level of the polarization signal in continuum 
windows is around 2-10"^ of the continuum intensity. 

We used a narrow-band G-band imaging channel for 
speckle-reconstructed context images and the TElecentric SOlar 
Spectrometer (TESOS: lKentischer et al.|[T998h for G-band spec- 



troscopy simultaneously with TIP. We aligned the speckle- 
reconstructed G-band images take n during the TI P maps with 
the same methods as described in Beck et al. l (l2Q07h :the TESOS 
spectra were not considered for the present study. 

Figure [T] shows overview maps of the observations: the 
aligned speckle-reconstructed G-band data, the continuum in- 
tensity Ic at around 1565 nm, the maximum linear polarization 
signal Lynax of the 1564.8 nm line, the maximum circular polar- 
ization signal Vmax^ the polarization degree /?, and finally, the po- 
larity of the magnetic fields that indicates if the field lines were 
(anti)parallel to the line of sight (LOS). The polarity was defined 
as ±1 (black and white in Fig. [T]) from the order of minimum 
and maximum Stokes V signal around the 1564.8 nm line for 
all locations with a polarization degree above the threshold dis- 
cussed below. Locations without significant polarization signal 
appear grey in the polarity map. Lmax^ ^max and p were defined 
as the largest polarization signal in a small wavelength range 
around the line core of 1564.8 nm in L = VQ^TTT^, Stokes 
y, and F - ^Q^ + + V^, respectively. Assuming that the 
speckle-reconstructed G-band data (leftmost image) is a close 
representation of the true solar surface structure, the TIP data 
show the same structure only worsened by the limitation of the 
spatial sampling, despite the 30 sec integration. The cadence of 
the G-band data was interrupted in both observations near the 
end of the scanning, and the G-band data has no exact timing 
information. It was thus impossible to obtain a good alignment 
to Ic for the last approximately IC of the scans. Lmax^ ^max and p 
are all displayed with an identical upper threshold of 1 % of Ic of 
polarization amplitude, p is dominated by the contribution from 
the circular polarization, i.e., nearly vertical magnetic fields, but 
there are several locations with strong linear polarization signal 
as well. Note that only a small fraction of the magnetic fields 
leads to bright points in the G-band (e.g. , x/y= 8-15/25 -20'' in 
the upper left image ), like also found bv ldeWiin et al.l 12008). 
The second map (lower row), taken around an hour later, shows 
again the strongest and stable magnetic flux concentrations of 
the previous observation displaced by around 5'' : the black patch 
in the polarity map at x/y = 38/55'' or the small white patch at 
18/58" can be re-identified in the 2nd map as well at 38/65" and 
18/63", respectively. 



3. Data analysis 

For a comparison with the study of LI08, we also calculated 
the total linear and circular polarization, Ltot and Vtot, by an 
integration of the absolute polarization signal L(X) and \V(X)\ 
over wavelength. These quantities are related to the longitudinal 
(Bl) and transversal (Bj) component of the magnetic field. The 
(anti)parallel direction of the longitudinal fields is introduced 
into Vtot again by a multiplication with the polarity map. Figure 
[21 displays maps of Ltot and Vtot- We added additionally maps of 
the ratio of linear and circular polarization, L/V. The ratio was 
calculated by an average over a few spectral pixels (-1 pm) 
near t he location of the largest V signal of 1564.8 nm (see lBeckl 
l2006h . The spatial pattern in these maps is strikingly similar to 
those found by LI08: the polarization signal shows "voids" of 
some granules extent (about 10" diameter in Ltot, 3-5"in Vtot) 
that are free from any clear signal; the linear polarization sig- 
nal appears in the shape of small "blobs" around the edges of 
the voids. The "blobs" o f enhanced L^^/^ are slightly smaller than 
individual granules (cf. Ilshikawa et al.1 ^2008). In the ratio LjV 
(right column), several more locations with stronger linear than 
circular polarization show up in the form of similar "blobs" as 
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Fig. 2. Left to right: Ltot, Vtot times polarity, ratio L/V. The red 
patches in L/V denote locations without any clear V signal to 
divide with. First and last map are displayed inverted. The blue 
rectangles in the upper left image denote two patches of large 
Ltot, between which in L/V (upper right image) another one 
shows up. Top/bottom: TIP Op. 001 and Op. 005. Scale is in arc- 
sec, grid lines have a 10 arcsec spacing. 



in Ltot- For example, at x/y=26/26" and 30/30" in the upper left 
image two "blobs" appear in Ltot (blue rectangles), whereas in 
L/V (upper right) a third one is additionally seen right between 
them at x/y=28/28". 

Figure [3] shows one example of slit- spectra from each ob- 
servation. There is a residual (vertical) fringe pattern in QUV 
of Op. 001 that comes from imperfect flatfielding , with a max- 
imum amplitude of around 0.0005 of Ic in Stokes U ; the am- 
plitude is smaller in Q and V. The amplitude of the fringes is 
actually so small that they do not show up even in the wave- 
length integrated map of Ltot (see Fig. [2]). The fringes are pre- 
sumably due to the glass window that protects the CCD chip 
inside the cryostat. We applied some additional corrections in 
the lower part of the slit from 0" to around 20", but point out 
th at it could still be im proved by making a PCA analysis like 
in iMartmez Gonzalez e t al. (2008c, a). For the 2nd observation, 
residual fringes are much weaker and can only be seen in Stokes 
Q. The inversion code is not able to reproduce the pattern due to 
the combination of its limited number of degrees of freedom, the 
least- square minimization it uses, and the extent of the fringes 
in wavelength. The characteristic wavelength of the fringe pat- 
tern is around two times as wide as one of the solar spectral 
lines; with the chosen setup, the inversion code has no degrees 
of freedom (see the next paragraph) to produce such a pattern. In 
this case, two synthetic profiles will show small deviations: the 
one that includes the fringe pattern (which the inversion code 
cannot generate), and the one that ignores the additional modu- 



Fig. 3. Example of top to bottom I QUV spectra for one scan step 
from Op. 001 at left and Op. 005 at right. The display threshold 
for QUV was set to ±0.001 h . 



lation. We thus assume that the influence of the fringes on the 
results is negligible. Comparing the two spectral lines at 1564.8 
nm and 1565.2 nm, it is clear that most of the linear polarization 
signals seen in the more sensi tive 1564.8 nm li ne are absent in 
the other, like already noted in Khomenko et all 12003): only the 
strongest linear polarization signals appear in both lines. Almost 
every profile along the slit shows a significant circular polariza- 
tion signal in both spectra. 

The noise level in the spectra is low enough that we could 
run a ''standard" inversion of the spectra with the SIR code 
(iRuiz Cobo & del Toro Iniestall 19921) . "Standard" means that we 
have not modified the inversion setup significantly to optimize 
it for the special type of data, to get a hopefully unbiased first 
guess of the magnetic field properties without enforcing a pri- 
ori knowledge. We used the same setup as B eck et al.l ('2007'): in 
each pixel we assumed one field-free and one magnetic atmo- 
sphere component, with additional stray light. We only modified 
the inversion to use three nodes in temperature for the magnetic 
component, and put equal weights for the polarization compo- 
nents 2, U, and V. The parameters of the magnetic field and the 
LOS velocities were taken to be constant with optical depth. The 
2-C+stray light setup mimics observations of unresolved mag- 
netic flux (mag. component) in equally unresolved intergranular 
lanes (field-free component), whereas the stray light corresponds 
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Fig. 4. Top: Example of inverted spectra of 1564.8 nm {left col- 
umn) and 1565.2 nm {right column). Black lines: observed spec- 
tra of top to bottom I QUV; red: best-fit profile of the inversion. 
The dash-dotted horizontal lines mark three times the rms noise 
level, the solid horizontal line the zero level. The vertical solid 
line denotes the rest wavelength. Bottom: polarization degree of 
1564.8 nm with the threshold for inversion {dash-dotted) and fi- 
nal rejection {horizontal solid line). 



to the granular contribution to the spectrcQ. Figures IH lA.ll and 
IA.2l show several examples of observed spectra together with the 
best-fit spectra retrieved by the inversion. We selected locations 
with a low polarization degree and mostly larger linear than cir- 
cular polarization signal for the figures. For all these examples, 
the signal in the less sensitive 1565.2 nm line already is close to 
or below the noise level, whereas the 1564.8 nm line still yields 
a "significant" signal that the inversion code can reproduce. 

For defining "significant", we used the polarization degree 
p. The routine that performed the inversion was using a thresh- 
old of 1 • 10"^ of Ic: if the maximal value of p of either 1564.8 
nm or 1565.2 nm was above this limit, the 2-C inversion with a 
magnetic component was done. As can be seen in Figs.l4lor lA.3l 
(right column, 3rd row), this level is reached already by (most 
probably) noise alone in two out of 13 cases for 1564.8 nm. We 
thus have set a slightly higher threshold for the analysis of the 
inversion result and rejected all profiles with p < 1.4 • 10"^ of 
Ic in 1564.8 nm. In both maps, around 63 % of the pixels are 



^ The average spectrum used as stray light contribution is blue- 
shifted and has unit continuum intensity. 



above the threshold that corresponds to 4 times the noise rms in 
/7( 3.4- IQ-^ of/.). 

ISheminoval (|2008) found that about 70 % of the surface area 
in a 2-D MHD simulation was covered with magnetic fields. 
OR07 give a fraction of 630 nm profiles of 87 % and 35.5 % 
for a polarization degree of 3 and 4.5 times the noise level of 
about 1 • 10"^ in their HINODE data. These numbers thus all 
agree that the area fraction covered by magnetic fields or show- 
ing polarization signal should be above 50 %, but Fig.[5]displays 
that the fraction of profiles with p > pthresh strongly depends 
on the threshold used. For comparison with a theoretical calcu- 
lation, we derived the same fraction of profiles with p > pthresh 
from spectra of 1564.8 nm synthe sized from a simulation run 
with the Co^bold code descr ibed inlSchaffenberger et al ] (l200l 
l2006l) and previously used in lSteiner et al.l (l2008l) . The snapshot 
was taken from the h20 run with an initial homogeneous hori- 
zontal magnetic field of 20 G. The simulation box corresponded 
to a 6! 6 square; its spatial resolution is much better than the 
present observations (cp. inlets 1 and 3). We thus have con- 
volved the simulations' spectra spatially with a Gaussian kernel 
of FWHM ~ a.'9. We also added a noise level of 2.2-10-4 to 
them. The original simulation spectra without convolution and 
noise {red dashed in Fig. [5]) have a much larger fraction of pro- 
files with a high polarization degree than the observations. After 
the convolution, the curve already matches roughly that of the 
observations {red solid). We ascribe the remaining diff'erence to 
the imperfect match of the spatial resolution, the diff'erent spatial 
sampling and the much larger solar area that enters in the obser- 
vations (40''x70'0. We only note that the convolved simulated 
spectra can be matched nearly perfectly to the observations by 
compressing the curve in the abscissa by a factor of 1.5 {blue 
dotted line). This compression can be understood as a reduction 
of the polarization degree that could be due to, e.g., the stray 
light that was not included in the convolution, or the usage of 
a Gaussian kernel with its steep drop with distance. The reduc- 
tion by 1.5 may sound quite large, but one has to consider that 
the eff'ect is to change, e.g., a 0.6 % polarization level to 0.4 % 
instead, for which only a small amount of additional stray light 
is needed. The slope of all curves changes clearly when reach- 
ing the noise level; the fraction of profiles is constantly close 
to 100 % for thresholds below the noise. We suggest that the 
fraction of profiles above a polarization threshold can be used 
easily to match the spatial resolution of simulation spectra and 
observations, and also is a good indicator of the spatial resolu- 
tion achieved in a specific observation. 




Fig. 5. Fraction of profiles with p > p thresh- Dash-dotted and 
solid vertical lines denote the inversion and final rejection 
threshold in the observations, respectively. The inlets show con- 
tinuum intensity maps of the original (1) and convolved simu- 
lation data (2), and a section of the same size from a TIP map 
(3). 
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Fig. 6. Comparison of the geometry in the calculation of SIR and 
a possible topology inside the resolution element, a, left to right. 
reference frame, LOS direction, and filling fraction, b: geometry 
in the calculation by SIR for field inclinations to the LOS of 
left to right 0, 25, and 90 degrees, c: possible solution for real 
topology. 



4. Inversion results 

The inversion code constructs its synthetic spectra, S , by the fol- 
lowing combination of inversion components: 

S = a-S stray + {I - o) ' {P • S magnetic + (1 - ft) ' ^ field-free) , (1) 

where S stray is the average profile of the whole map that is taken 
to describe the stray light contribution and a its fractional con- 
tribution to the spectrum. /3 and 1 - are the relative contribu- 
tions of the magnetic and the field-free component. Assuming 
that S stray simply corresponds to another field-free component 
with a fixed velocity, we define the filling fraction / of the mag- 
netic fields by / = {I - a) • JS. The additional quantities from 
the inversion code are the field strength, B, field inclination to 
the LOS, 7, field azimuth, i//, LOS velocities for each compo- 
nent, and the temperature stratification. The observations were 
on disc centre, thus the LOS inclination y is identical to the field 
inclination to the surface normal. 



From the inversion results, we calculated an estimate of mag- 
netic flux by 

= /.5.A, (2) 

where A is the area corresponding to a single pixel of 0'.'5 x 0^.^35. 

For this estimate of the magnetic flux, we did not consider 
the field inclination. We motivate this with the assumption that 
the volume contributing to the spectra from one pixel corre- 
sponds roughly to a (275 km)^ cube. The spatial extent in the 
horizontal dimensions {dx, dy) is 360 km by 250 km. For the 
extension of t he formatio n height, H, of the IR line at 1564.8 
nm, Cabrera Solana et al. 1(12005) sive a range of around two dex 
in log T, e.g., from around to -2 for velocity perturbations or 
from 0. 5 to -2.5 for magnetic fields. In the HSRA model atmo- 
sphere (iGingerich et al.lll97ll) they used, this range corresponds 
to around 250 km. Figure [6] shows one possible solution how 
the calculation of SIR, which determines a filling fraction in the 
horizontal plane and then solves the radiative transfer along the 
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Fig. 7. Inversion results of Op. 001 (top row) and Op. 005 (bot- 
tom row). Left to right: field strength, LOS inclination, filling 
fraction, total flux. Red contours in O trace strong linear polar- 
ization signals. The display ranges are given by the labels of the 
color bar at bottom. Pixels not inverted are shown black. 



LOS, could be related to the real topology inside the pixel. We 
assumed a constant filling fraction of 25 % of the horizontal area 
Axy and plotted three cases for field inclinations to the LOS of 
0, 25, and 90 degrees. To obtain the magnetic flux, the filling 
fraction has to be multiplied with the corresponding area per- 
pendicular to the field direction. For horizontal magnetic fields, 
the respective area is, however, not Axy but Axz or Ay^, depending 
on the field azimuth. If dx - dy - H, the surfaces of the cube are 
equal and thus the field orientation can be ignored to first order. 
This is only approximative, as for example in the 90-deg-case, 
SIR calculates the synthetic profile assuming that the horizontal 
fields are present everywhere inside the formation height, which 
not necessarily will give the same result as if they fill only half 
the formation height with doubled filling fraction. 

Figure [7] shows maps of B, y, / and O for the two observa- 
tions. Magnetic fields with kG or more are restricted to vertical 
fields (7^^ and 2nd column, y ~ or 180°). The kG fields show 
a high magnetic filling fraction, but also parts with horizontal 
fields can reach the same level in /. The more inclined fields, 
however, do not appear prominently in the flux map (rightmost 
column). We overplotted contour lines tracing the strongest lin- 
ear polarization signal (cf. Fig. O onto the flux map. Locations 
of large magnetic flux are almost exclusively related to vertical 
fields with little linear polarization signal. 

The total amount of magnetic flux in both maps calculated 
with Eq. ^ is around 6-10^^ Mx. There is a slight imbalance 
of 53 to 47 % for magnetic flux parallel or antiparallel to the 
LOS. Including the field inclination in the calculation to obtain 
the LOS magnetic flux by multiplying O with cos y, the ratio 
is reversed to 47 to 53 %. As the fraction is close to 50 %, we 
think that the FOV gives a sufliiciently large statistical sample of 
quiet Sun regions where no strong imbalance in favor of a special 
field orientation is to be expected. The histogram of inclinations 
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Fig. 8. Histograms of (clockwise, starting left top) magnetic flux, 
filling fraction, field inclination, field strength. Blue and red lines 
correspond to the least/most inclined fields, black lines to the full 
FOV. 



using the full range from to 180 deg was also fairly symmetric, 
so we have in the following converted field inclinations above 
90 degrees to their corresponding values in the range from to 
90 degrees. 

Figure [5] shows histograms of B, O, y and /. We calculated 
the histograms for three diff'erent samples there: the full FOV 
(black), all locations with nearly vertical magnetic fields (blue, 
y < 25°), and nearly horizontal magnetic fields (red, y > 75°). 
For the histogram of y we separated into magnetic fields above 
and below 500 G instead. The visual impression of the 2-D maps 
is confirmed by the histograms: the nearly vertical fields corre- 
spond to the locations with large magnetic fluxes and high field 
strength (300-1000 G); the nearly horizontal fields are limited to 
below around 500 G. 

The distribution of magnetic field strength for the full 
FO V agrees wit h previous observations of the IR lines 
( Kho menko et al.i r2003: Martinez Gonzalez et al. 2008c) or the 
field strength distribution derived by Sheminova ( 2008) from 2- 
D MHD simulations, especially in the point of the most probable 
field strength of around 250 G. We point out that a) the "bump" 
at around 1.5 kG that commonly is found from the analysis of 
630 nm data (e.g., Dominguez Cerdena et al. 2006a, OR07) is 
missing, and that b) the decrease of relative occurrence for fields 
below 250 G is not due to the polarization detection limit. Figure 
IA.2I shows 3 examples of weak fields below 200 G (2nd row) 
that still could be reliably reproduced by the inversion. The spa- 
tial resolution of observations is also crucial for the detection 
of these weak fields (see Fig. [5]). The distribution of the nearly 
horizontal fields (red curve at bottom left of Fig. [8]) shows an 
enhanced fraction of these weak fields below 200 G which im- 
plies that the detection limit due to the noise is not yet reached 
at 200 G. In the question of the seer ningly preference for kG 
fields in the analysis of 630 nm data, iMartinez Gonzalez et al.l 
(l2006bl) have cautioned that the thermodynamics have to be 
treate d consistently for t he 630 nm line pair in the weak field 
limit. iRezaei et al.l (1200 7) were able to derive a distribution with- 
out the bump using a SIR inversion of 630 nm data of lower 
spatial resolution than the nowadays available HINODE obser- 
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Fig. 9. The same histograms as in Fig.[8]but separated into gran- 
ules {black lines) and intergranular lanes {red). The inlets show 
the masks for the separation. 



vations. It would be worthwhile to analyze the HINODE data 
without employing the Milne-Eddington (ME) approximation to 
see if the field strength distribution changes with a more sophis- 
ticated treatment of the thermodynamics. 

The field inclination shows some clear trends: for fields be- 
low 500 G, an inclination above 45° is preferred with a plateau of 
constant probability up to 90°, whereas strong fields above 500 G 
are restricted to an inclination of 35° or less. The distribution for 
the full FOV shows a small peak near 0° inclination (= the very 
small area fraction of "network" fields inside the FOV), and then 
a rapid increase of occurrence up to 15° followed by a slow re- 
duction towards 90°. The shape of the distribution matches more 
closely to the on e given by QR07 in their Fig. 3 for the MHD 
simulation run of lVogler etal.1 (120051) than the one these authors 
derived from their observations. 

A second method of classification into diff'erent type of struc- 
tures is to separate the FOV into granules and intergranular lanes 
(IGLs). We derived a mask using the continuum intensity maps 
and a threshold of 0.995 and 1 .005 of Ic to define IGLs {black in 
the inlets in Fig. [9]) and granules {white), respectively. Locations 
with an intensity between the two values were excluded {green). 
The masks retain the granulation pattern of the continuum in- 
tensity. At the near-infrared wavelength of the observations and 
at their spatial sampling, magnetic elements do not show up as 
bright in the continuum on disc center, and thus should not in- 
cidentally have been counted as part of the granular sample. We 
then calculated the histograms of the same quantities as in Fig.[8l 
The pattern is very similar to that in Fig. [8l in the IGLs, one 
has a predominance of more vertical, stronger fields with larger 
magnetic flux. A strong difl'erence between granules and IGLs is 
seen in the magnetic filling fraction. IGLs show a larger contri- 
bution of high filling fractions centered at 0.2-0.3, whereas the 
granules show a distribution that peaks close to / = and drops 
steeply for larger /. This indicates that the magnetic flux in (or 
above) granules is more difl'use than in the IGLs and gets rather 
dispersed than concentrated, whereas in the IGLs the converging 
granular flows advect magnetic flux and concentrate it. 

The relative amount of magnetic flux in concentrated ver- 
tical fields in IGLs and the diff'use weak horizontal fields in 




Fig. 10. Scatterplots of magnetic flux vs inclination {top) and 
field strength vs inclination {middle). Red crosses are binned 
data. Bottom: total magnetic flux in the full FOV as function of 
field inclination in 1 deg bins. The black vertical line marks 45°; 
the numbers give the total magnetic flux of all fields below or 
above this limit. 

granules is of strong interest, e.g., for reconciling Hanle and 
Zeeman measurements with low spatial resolution that yielded 
quite diff'erent average magnetic flux densities, with Zeeman- 
based observations lower by around an order of magnitude. The 
histograms of Figs. [8] and [9] indicated a tight relation between 
field strength, magnetic flux and field inclination. We have thus 
also created scatterplots of inclination vs the other two quantities 
{top two panels of Fig. [TO]). We remark that the field inclination 
did not enter into the calculation of the magnetic flux, but only 
filling fraction and field strength. In the scatterplots, there is a 
clear trend that the inclined fields above 45° are limited to be- 
low 500 G and below around 5-10^^ Mx. Larger fluxes are only 
achieved by nearly vertical fields (cp. Ishikawa & Tsuneta 2009^, 
their Fig. 6). To obtain the total amount of magnetic flux as func- 
tion of inclination, we added the magnetic fluxes of all locations 
with a field inclination in [7,7+ 1°] {lower panel). The distribu- 
tion is similar to the scatterplot at top; we find that around 2/3 of 
the total magnetic flux comes from magnetic fields with 7 < 45°. 

In observations with low spatial resolution, the weaker and 
more inclined magnetic fields could not be detected because 
their linear polarization signal was spatially smeared out to val- 
ues below th e noise level (see Fig. [TTl o bservations are described 
in det ail in iKhomenko etaD|2Q05i iMartinez Gonzalez et al.l 
l2008ch . These data of 2003 were taken in the same spectral range 
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Fig. 11. Linear polarization map Ltot of an observation of quiet 
Sun with TIP @ 1.56 micron from 17.8.2003. Tick marks are in 
arcsec. 



here and at 0.2 in OR07. The field inclinations diff'er most: we do 
not find an increase in the distribution for horizontal fields, but 
a probability declining slowly with inclination for fields with 
inclinations above 20°. OR07 and LI 08 both found a clear dom- 
inance of horizontal over vertical fields and transversal over lon- 
gitudinal flux, respectively. With the same caveat than in LI08 
that the vertical and horizontal extent of the transverse fields ac- 
tually cannot be derived at once from the spectra, we found that 
2/3 of the total flux comes from fields with inclinations below 
45°, whereas LI08 give a ratio of 5 in favor of the transversal 
fields. 



as the present investigation, with an identical integration time of 
30 sec per spectrum, but an image correction with only a corre- 
lation tracker to remove image shifts during the exposure. A sin- 
gle position at (48'',22'0 in a 60''x35'' FOV shows significant 
linear polarization signal that would indicate inclined magnetic 
fields. iKhomenko et al.l (l2003l) . however, were able to detect lin- 
ear polarization signals in 50- sec integrated data from TIP taken 
on July 29th 2000 and September 5th 2OOC0. They found that 
around 20 % of the pixels inside the FOV showed a significant 
linear polarization signal that could be analyzed quantitatively. 



Reliability of proxies in magnetic flux determination As we 

were able to successfully run an inversion on most of the pix- 
els due to the low noise level of our data, we have the informa- 
tion on field strength, field inclination and magnetic flux with 
some certainty. This allows to test the validity of some of the 
approximations in LI08 who convert the total linear and circular 
polarization signal into their estimates of transversal and longi- 
tudinal flux. To this extent, we calculated the longitudinal (O^) 
and transversal (Oj^) flux from the inversion results by 



= / • 5 • cos 7 • A 
Oj^ = / • 5 • siny • A . 



(3) 



5. Discussion 

General properties of the QS magnetic fields The weak 
horizontal magnetic field of the solar internetwork regions 
is only accessible to observations of high spatial resolu- 
tion and high polarimetric sensitivity. The presence of these 
fields has been expected from Hanle measurements, and 
traces o f them have been found in some observations before 
(iLites et al. 1996; Khomenko et all |2003|; 'Harvey et al. 2007; 
iMartinez Gonzalez et al. 2008b; c). They have only recently been 
seen prominently with dat a from the spectro polarimeter on- 
board the HINODE satellite (iKosugi et al.ll2007h . as reported by 
OR07 and LI08. From a comparison of the present Fig. [2] with 
Figs. 2 and 7 of LI08 we conclude that the present observations 
trace the same solar structures, i.e., ground-based observations 
aided by real-time adaptive optics correction can reach the nec- 
essary spatial resolution during spells of excellent seeing condi- 
tions, and hence, the results of the analysis should be compatible. 

Some of the present findings are at odds to the HINODE 
data, but we think that we can trace back most of these con- 
tradictions to the analysis methods, the diff'erent magnetic sen- 
sitivity of the spectral lines, and finally, to a physical reason. 
The distributions of field strength show again the mismatch be- 
tween IR and VIS found in previous studies: a significant frac- 
tion of kG fields in VIS (OR07) that do not appear for IR 



As first order approximation, the circular polarized Stokes V 
signal should scale with Ol, whereas the linear polarization 
L = VG^TTT^ should scale with (jJelferies et al.lll989l) . 
From the observed profiles, we derived the maximal linear and 
circular polarization of the 1564.8 nm line {Lmax and Vmax, see 
Sect. O and the total linear and circular polarization {Ltot and 
Vtot, see Sect. O. In the derivation of Ltot, the noise in both Q 
and U enters. Following LI08, this contribution can be mini- 
mized by rotating the linear polarization signal to the "preferred 
azimuth frame" where the signal is concentrated in Stokes Q (see 
Appendix |B]). We will thus in the following use 2^^^ (rotated) = 
J |2(/l)|J/l after the rotation of the spectra as a measure of the 
linear polarization instead of Ltot- Figure [121 shows scatterplots 
of Or vs VGmax and ^[Q^t, and vs Vmax and Vtot- We also 
added a plot of field strength vs ^|Qtot and Vtot- It can be seen 
that there is on average a clear relation between (^t or and the 
corresponding quantities derived directly from the profiles {red 
points), but the scatter around this relation is large (around 30 % 
diff'erence between mean and minimum or maximum values). 

The polarization amplitudes, and hence, also the integrated 
polarization amplitudes measure the ratio of polarized to un- 
polarized light. Their absolute values strongly depend on the 
thermodynamical structure of the atmosphere and the resolution 
of structures. This gets even worse for VIS lines in the weak- 
field limit where the field strength also changes the amplitudes 
of polarization lobes, not only their location in wavelength. In 
Appendix |B] we tested the derivation of calibration curves be- 
tween transversal magnetic flux and Qtot under varying assump- 
tions (constant inclination, constant flux, temperature variation). 
We find that only varying the temperature stratification at con- 
stant magnetic flux is almost sufficient for producing the scatter 
seen in Fig. [121 In total, the scatter of the relationship and the 
different calibration curves that can be obtained by, e.g., only 
varying the inclination by 10 degrees put the derivation of mag- 
netic ffux from proxies as Qtot or Vtot into some doubt. 

However, besides from the observations with HINODE that 
favor a dominant fraction of horizontal fields recently also 
some simulation results have suggested more hor i zontal mag- 
See http://www3. kis^m^freiburgjd^ then previously. ISchiissler & Vogle3 (l2008l) and 



data. iBellot Rubio & Conadosl j2003l) . iMartmez Gonzalez et al.l 
(I2006bl) . and Rezaei et al. (2007b have shown that this feature is 
related to the noise level of the data and the thermodynamics, 
because IR and VIS spectra to gether can be reproduced sim ulta- 
neously with weak fields ( Mar tinez Gonzalez et al.ll2b08cl) . The 
usage of a ME method in OR07 possibly has also a contribution 
to this effect. The fie ld st reng th distr ibution here peaks at around 
250 G, whereas ( Orozco Suarez et al.|[2007l) give 90 G. This can 
again be related to the weak-field limit of the 630 nm lines to- 
gether with the ME approximation that makes the VIS lines less 
reliable than the stronger splitting IR lines. The magnetic filling 
fraction / compares fairly well to one minus the stray light con- 
tribution of OR07; magnetic filling fractions peak at around 0.1 
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Fig. 12. Scatterplots of transversal and longitudinal 
flux vs maximal linear (left top) and circular po- 
larization degree (left bottom), same vs total linear 
(middle top) and circular polarization (middle bot- 
tom), and scatterplot of field strength vs total lin- 
ear (right top) and circular polarization (right bot- 
tom). Black: all data points; red: after binning. The 
blue line in the middle top panel is a linear fit to the 
binned data. 



ISteiner et"an (12008'. ST08) have reported the ratio of horizontal 
to vertical fields in their respective simulations and give numbers 
from 2 to 6, in accordance with the factor of 5 derived by LI08. 
This seems incompatible with the present investigation, but both 
simulations show that the ratio is strongly dependent on height 
in the atmosphere. ST08 show in their Fig. 1 the field strengths 
as function of height for two simulation runs (vlO=initial 10 G 
vertical field, h20= 20 G horizontal). At a height of around 200 
km, which is appropriate for the IR lines at 1565 nm used here, 
the ratio of vlO is actually in favor of the longitudinal fields, 
Bhor/Byer ~ 0.82. This is significantly closer to the ratio of 
0^>45o /0^<45o = 0.42 than a factor of 5 that ST08 derive after 
synthesizing 630 nm spectra of their simulation and analyzing 
them in the same way as LI08. 

Vertical unsigned flux density Various authors have de- 
rived the average unsigned flux density of the verti- 
cal mag netic fields from Zeeman sensitive spectral lines 
(iKeller e t al. 1994; Lites 2002; Dominguez Cerdena et al. 2003; 
'Khomenko et al. 2005a; Dominguez Cerdena et al. 2006a, b; 
N Orozco Su a rez e t al. 2007; Martinez Gonzalez et al. 2008c; 
Lites et al.1 [20081: ICarroU & Ko pf 2008). In the older literature, 
magnetic field proxies like the Stokes V amplitude or integrated 
area were used; in the most recent ones an inversion of spec- 
tra. The obtained values for {3^} scatter strongly, depending 
on the observations used or the analysis method employed, but 
in generally yield values between 5 to 30 G. LI08 in addition 
give a value of around 50 G for the transversal magnetic flux. 
Observations using the Hanl e eflTect yielded sign ificantly larger 
values of around 100 G re.g.. lTrujillo Bueno et al.ll2004D . 

Figure [13] shows the distribution of {B^) in our observations, 
derived from the longitudinal flux by division with the area A 
(cf. Eq. ([2])); the bin size was 0.5 G. The distribution shows a 
maximum near 2 G, but otherwise is in good agreement with 
a monotonically increasing probability for weaker fluxes (solid 
red line, exponential decay law). The diflTerence between the ex- 
ponential decay and the observed distribution for O < 2 G ac- 
tually only covers a small area of far less than 10 % relative 
frequency (red filled area in the upper panel of Fig. [13] from 
to 2 G, see the inlet). In the double-logarithmic display in 
the lower half, we displayed the first 7 histogram points with 
crosses. It can be clearly seen that from high fluxes up to the 
maximum the points follow well the exponential law; only the 
first four points referring to the weakest fluxes below 2 G devi- 
ate. Taking into account that the noise imposes a hard thresh- 



old for the detection of the smallest magnetic fluxes, we believe 
that the occurrence of the maximum is fully artificial. Very weak 
fluxes produce smaller and smaller polarization signal that can 
escape detection c ompletely or are eventually reje cted by the in- 
version threshold. iPietarila Graham et al.l (l2009l) claim that the 
distribution of the magnetic flux density cannot be determined 
from observations because of the peaked appearance of the prob- 
ability distributions like in Fig. [13] However, given the actually 
small deviation between the observed distribution and an expo- 
nential decay law predicted by MHD simulations, we think that 
the probability distribution function of magnetic flux can be de- 
rived from observations when one takes into account that the 
weakest fluxes must escape detection due to the noise level. 

Using the strongly Zeeman- sensitive IR lines at 1.56 jim, it 
is possible to disentangle the field strength B on the part of the 
surface that actually is covered by fields from the area filling 
fraction / of the almost certainly unresolved magnetic struc- 
ture. This makes an important diflTerence for the calculation of 
the magnetic flux density, compared with the spectral lines at 
630 nm used by the HINQDE spectropolarimeter and in the 
work of IPietarila Graham et al.l (120091) . The 630 nm lines are in 
the weak field limit for the weak magnetic fields of the quiet 
Sun, i.e., they suffer from a strong interplay between thermody- 
namics and magnetic field properties, as also demonstrated in 
Appendix |B] It is more difl&cult to determine B and / accurately 
for these lines, even if one takes care about the thermodynamics 
(iMartinez Gonzalez et alll2006lj) . 

Normalizing the observed magnetic flux to the complete 
FOV, including the locations with no significant polarization sig- 
nal, the average values for the total, longitudinal and transversal 
flux using Eqs. ([2]) and ([3]) are 22 G, 16 G, and 1 1 G, respectively. 
If one considers only the area with actually observed polariza- 
tion signal, the corresponding values are 26 G, 20 G, and 13 G. 
The values are on the upper range of Zeeman base d measure- 
ments. As mentioned by, e.g.. IPietarila Graham et al.l (120091 the 
cancellation of Stokes V signals due to unresolved opposite po- 
larities inside a resolution element causes underestimation of the 
real amount of magnetic flux. Our observations with a O'.' 5x0^.^37 
sampling are prone to this. The average flux densities thus can 
only be taken as lower limit for the total magnetic flux, but we 
believe that they actually give a rather solid limit for the total 
magnetic flux observed at the V resolution achieved. 
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Fig. 13. Histograms of vertical flux density in linear (top) and 
logarithmic scale (bottom). The vertical lines mark 2 G. Black 
and red lines show the observations and an exponential decay 
law. The inlet in the upper panel shows a magnification of the 
low flux range. 



6. Summary & conclusions 

We have analyzed two long-integrated data sets of quiet Sun 
regions obtained under extremely good seeing conditions with 
the TIP at 1.56 jim. We find that the data allow to study the 
weak horizontal internetwork fields that appeared prominently 
in HINODE observations recently (OR07 and LI08). Thanks to 
a lower noise level than in the HINODE data and the higher 
sensitivity of the IR lines, we could run an inversion of the 
maps with SIR that saves us from using strong approximations 
in the derivation of magnetic field properties. We find that the 
more inclined/"horizontal" fields (y > 45°) are generally weaker 
(<500 G) and contain in total less magnetic flux than the less 
inclined/'Vertical" fields. As the data properties (spatial scales 
of observed structures, pixels with significant polarization sig- 
nal) and the inversion results (magnetic filling fraction, relation 
to granules/IGLs) seem roughly to match to the HINODE data 
and their analysis in OR07 and LI08, we ascribe the rever- 
sal of the ratio of transversal to longitudinal fields (HINODE 
data: 5, IR data: 0.42) to a) the diff'erent formation height of the 
spectral lines and b) the analysis methods (ME inversion, cal- 
ibration from integrated polarization signal to magnetic flux). 
We conclude that the issue should best be investigated with si- 
multaneous IR/VIS data, as simulations pre dict a strong depen- 
dence of the ratio on geometrical height ( Steiner et al I l2008l) . 
Whereas the existence of a significant fraction of inclined mag- 
netic fields is confirmed by the present IR observations, their 
influence on the dynamics of the solar atmosphere is not clear 
yet. They do not seem to have an impact on the granulation pat- 



te rn, but could contribu t e to ch romospheric heating as discussed 
in llshikawa & Tsunetal (l20Q9h . This issue can be addressed by 
a combination of repeated spectropolarimetry of photospheric 
fields with chromospheric diagnostics. 
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Appendix A: Examples of inverted spectra 

Figures lA.ll and IA.2I show several profiles taken from the first 
and second long-integrated observation on 21/05/08 (Op. 001 
and Op. 005), respectively. The positions of the profiles are 
marked with a consecutive number in Fig.[T] The profiles shown 
were selected to have a small polarization degree that in some 
cases was barely sufficient to meet the inversion threshold (e.g., 
profiles no. 6 and 7). Below the spectra, the temperature strati- 
fications that were used for the generation of the best-fit spectra 
are shown. With 3 nodes in temperature, the SIR code can use a 
parabola for changing the stratification; the parabola shape ap- 
pears quite prominent for many of the locations. Note, however, 
that the IR lines at 1 .56yum are not s ensitive to the temperature i n 

the atmosphere above log r 1.5 (ICabrera Solana et al ]l2005l) . 

Only one profile corresponds to a kG field (Fig. lA.2l top middle, 
no. 8). Figure |A3] shows the polarization degree of 1564.8 nm 
for all profiles of the previous figures. Profile no. 9 exceeds the 
inversion threshold of 0.001 of near -1-750 mA with a spike 
that presumably is not of solar origin, but noise in the Stokes U 
profile. The final rejection threshold of 0.0014 is, however, only 
reached by signals clearly related to the Zeeman eff'ect (multiple 
double or triple lobes). 

Reliability of the inversion results As discussed in Sect. 3, we 
used a constant value for the magnetic field strength (5), inclina- 
tion (y), azimuth {ij/), and the LOS velocity. This inversion setup 
cannot reproduce antisymmetric Stokes 2 or or symmetric 
Stokes V profiles that would require gradients in the magnetic 
field strength and the velocities along the LOS. The inversion 
was initialized with the same model atmosphere on all pixels 
{B = 0.9 kG, i// = 65 deg), only the inclination was modified 
to be 10 or 170 deg depending on the polarity. In the inversion 
process, the equal weight used for QUV in the calculation of 
naturally favors the component with larger polarization signal. 
For example, in profile no. 5 in Fig. lA. li the Q and U signals are 
larger than the V signal by almost an order of magnitude, lead- 
ing to a better fit quality for Q and U than for V. In polarimetric 
data of lower S/N, a diff'erence of this order usually implies that 
the weaker signal is not seen at all. 

SIR calculates an error estimate for the free fit parameters 
using the diagonal elements of the covariance matrix, approxi- 
mated by response functions (B ellot Rubio et al. 2000, SIR man- 
ual). The error estimate depends on the number of degrees of 
freedom in each variable; for parameters constant with optical 
depth thus a single value is returned. The error estimate gives, 
however, only the information on how reliable the best-fit so- 
lution of the found -minimum is inside the chosen inversion 
setup. The estimated errors for the profiles shown in Figs. 4, 
lA.ll and IA.2l are noted on the Stokes / panel of Fe 1 1565.2 nm 
line. The average uncertainties in the calculated magnetic field 
strength and inclination angle given by SIR are ±50 G and ±10 
d eg, respective ly. The values agree with a previous error estimate 
in lBeckl(l2006h that was derived from a direct analysis of the pro- 
file shape of the 1.56 jim lines (Table 3.2 on p. 47; SB 50 G 
and ~ 5 deg). 



the "preferred azimuth frame" where the linear polarization sig- 
nal is concentrated in Stokes Q. To achieve this, we determined 
the azimuth angle from the ratio of U and Q and rotated the 
spectra correspondingly to maximize the Stokes Q signal. The 
scatterplot in Fig. IB. 1! compares the previously used total linear 
polarization, Ltot = J + U^(A)dA, with the corresponding 
Qtoti^otSited) = J |2(/l)|J/l as a measure of the linear polariza- 
tion. The rotation of the spectra reduces the noise contribution 
by a constant amount, but the old and new values otherwise have 
a linear relationship with a slope close to unity. 

We then averaged the rotated Q spectra over all spatial posi- 
tions exceeding the polarization threshold for the inversion. The 
average Stokes Q spectrum is used as a spectral mask by LI08, 
but unfortunately their method fails for the infrared lines. The 
wavelengths around the line core have negative values in the av- 
erage Q profile (Fig. IB. 2!) which prevents using it in the same 
way as in LI08. We thus used 2^^^ (rotated) as defined above in- 
stead which we think to be equivalent to that of LI08 despite of 
not using a (somewhat arbitrary) spectral mask. 

The plot of 2^0^ (rotated) vs transversal flux (Fig. [121 middle 
upper panel), showed considerable scatter that already put the 
use of a single calibration curve into doubt. We thus not only 
tried to obtain a calibration curve, but also to quantify the efl'ect 
of various parameters on the obtained relation. The upper part 
of Fig. |B3l shows calibration curves of Qtot vs field strength for 
diff'erent field inclinations y. The uppermost curve for y = 90° 
corresponds to the one used by LI08. With the assumption that 
the field inclination not necessary equals 90°, one already finds 
that one and the same value of Qtot can be obtained in a range 
of around 200-550 G in field strength. The same eff'ect is shown 
in the middle part where the magnetic flux, O = 5 sin y, was 
kept constant at 1.8-10^^ Mx, B was varied, and y was derived 
accordingly from y = arcsin(0/5). Again a range of around 200- 
500 G in 5 gives the same value of Qtot- As final test we chose 
to investigate the influence of the temperature stratification on 
the resulting Q^^rvalue. We kept magnetic flux, field strength 
and field inclination constant at (1.8-10^^ Mx, 20 G, 75°), and 
synthesized spectra for diff'erent temperature stratifications. We 
used 10000 temperature stratifications that were derived for the 
magnetic component in the inversion, and thus can be taken to 
be an estimate for the range of temperatures to be expected in 
the quiet Sun. The histogram of the resulting 2^^^ -values is dis- 
played in the bottom part of Fig. IB. 31 ^/Qtot ranges from nearly 
zero up to 0.01, which roughly also corresponds to the scatter of 
^JQtot in Fig. [121 We thus think that the biggest contribution for 
the scatter comes from temperature eff'ects. We remark that we 
used a magnetic filling factor of unity in all calculations. Any 
additional variation of the filling factor for unresolved magnetic 
structures will increase the scatter of Qtot even more. 

We conclude that the usage of a calibration curve for a 
derivation of transversal magnetic flux from Ltot or Qtot, regard- 
less of the exact calculation of the wavelength integrated quanti- 
ties, is not reliable for a solid estimate, mainly because the strong 
influence of the thermodynamical state of the atmosphere on the 
weak polarization signals. 



Appendix B: Calibration of Ltot to transversal flux 

We tried to follow the procedure described in LI08 to obtain 
a calibration of linear polarization signal into transversal mag- 
netic flux that is independent of the inversion results. To re- 
duce the influence of noise, the latter authors first determine 
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Fig. A.l. Examples of observed (black line) and best-fit spectra (red) in Op. 001 ( 1st and 3rd row). The dash-dotted horizontal 
fines in QUV mark three times the rms noise level, the solid horizontal line the zero level. The vertical solid line denotes the rest 
wavelength. The 2nd and 4th row show the corresponding temperature stratifications of the magnetic component (solid), the field- 
free component (dash- dotted), and the HSRA atmosphere that is used as initial model (dashed). Field strength and LOS inclination 
and their respective errors are given in the plot of Stokes / of 1564.8 nm (upper left in each panel); the number of each profile is 
given in the upper left corner of each panel. 




Fig. A.2. Same as Fig. [Allfor TIP Op. 005. 
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Fig.A.3. Polarization degree of 1564.8 nm for the profiles shown in Fig. \A.l\ (upper two rows) and IA.2I (lower two rows). The 
dash-dotted and solid horizontal lines denote the inversion and final rejection threshold, respectively. 




Fig. B.l. Scatterplot of the integrated Stokes Q signal in the pre- 
ferred reference frame vs the total linear polarization without 
rotation. Solid line: unity slope; dashed line: unity slope with an 
offset of 0.0001. 
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Fig. B.2. The average Stokes Q profile. 
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Fig.B.3. Top: calibration curves from Qtot into field strength B 
for field inclinations y from 10 to 90 deg (bottom to top). The 
horizontal dotted line is at Qtot = 0.007; the solid part of it de- 
notes a range in B that gives the same Qtot at diff'erent y. Middle: 
^JQtot vs field strength for constant magnetic flux. Dotted line 
and solid part as above for ^JQtot = 0.0063. Bottom: histogram 
of ^/Qtot for constant flux but varying temperature stratifications 
T. The vertical line denotes the value resulting from the HSRA 
atmosphere model. 



